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Abstract 

One of the remaining important problems of life cycle inven¬ 
tory analysis is the allocation problem. A proper solution of this 
problem calls for a proper understanding of the nature of the 
problem itself. This paper argues that the established definition 
of the allocation problem as the fact that one unit process pro¬ 
duces more than one function, is not appropriate. That defini¬ 
tion points to an important reason of the occurrence of the prob¬ 
lem, bur the situation of internal (closed-loop) recycling already 
indicates that there may be product systems which contain multi¬ 
function processes, but which nevertheless need not exhibit an 
allocation problem. The paper proceeds by examining a number 
of simple hypothetical cases, and proposes a precise and opera¬ 
tional definition of the allocation problem. This enables a sys¬ 
tematic categorization of approaches for dealing with the allo¬ 
cation problem. 


Keywords: Allocation problem, algorithmetic procedures; LCA 
methodology, allocation problem, algorithmetic procedures; ma¬ 
trix inversion, allocation problem, LCA methodology 


1 Introduction 

Many books and papers which discuss the principles of life 
cycle inventory analysis describe the allocation problem as 
the fact that one process provides more than one valuable 
function (see, for instance, Consoli et al., 1993; Huppes & 
Schneider, 1994; Lindfors et al., 1995; Anonymous, 1998 and 
Azapagic, 1996). This then immediately offers a starting point 
to discuss solutions to the allocation problem: the multi-func¬ 
tion process must be split into a number of virtual processes 
by means of so-called allocation factors, or a process must be 
subtracted from the product system so that the coproduced 
function avoids that some other system produces it. 

This way of defining the allocation problem has proven to 
be an efficient way of communicating: the allocation prob¬ 


lem is now perceived by many scientists and practitioners in 
more or less one meaning. Also, solutions to the allocation 
problem have been brought to a much higher level than be¬ 
fore. The earlier books wrote about a simple allocation on a 
mass basis (e.g. Faya et al., 1991; Fecker, 1992), while the 
current developments show a sophisticated hierarchy of so¬ 
lutions, sometimes with criteria for accepting or rejecting a 
certain solution (e.g. Huppes Sc Frischknecht, 1995; Lindfors 
et al. 1995; Anonymous, 1998). 

Nonetheless, there remain some unclarities in these discus¬ 
sions. One instance of this is provided by the situation of 
internal recycling (alternatively called closed-loop recycling): 
here we definitely have a multi-function process, providing 
the functions waste treatment and production of secondary 
material, but there is not necessarily an allocation problem. 

As an example, consider the simple system of the two proc¬ 
esses "electricity production” and "fuel production", with 
the electricity production process producing a waste that 
can be used as an input for the fuel production process. The 
fuel production process thus effectively recycles waste into 
fuel. The complete specifications are: 

• electricity production: output of 10 kWh electricity, out¬ 
put of 2 kg waste, input of 1 1 fuel; 

• fuel production: output of 100 1 fuel, input of 500 kWh 
electricity, input of 200 kg waste. 

As a functional unit we choose 1000 kWh electricity. It is 
easily checked that a multiplication of the electricity pro¬ 
duction process by 200 and a multiplication of the fuel pro¬ 
duction process by 2 leads to a consistent balance: there is 
a gross production of 2000 kWh electricity, half of which is 
used as an input for fuel production. The 200 1 fuel that is 
produced is exactly consumed in the production of electric¬ 
ity, and exactly the 400 kg waste that is produced is fed 
into the production of fuel. In other words, we have in¬ 
cluded a recycling process without having performed an 
allocation step. 
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From the above example, we can easily illustrate an impor¬ 
tant conclusion of the mathematical portion of this paper. A 
central task of the inventory analysis is the scaling of the 
data of the unit processes such that there is a system-wide 
balance per material and per service. If there are / of these 
balances and J unit processes, we need to find J scaling fac¬ 
tors by means of / equations. Finding the scaling factors is 
straightforward if / = /, that is, if all processes produce a 
single function. When some processes produce more than one 
function, we have more equations than unknowns (I > J). This 
overdetermined system of equations will in general not have 
a unique solution, but it may have one in special cases, among 
which is the case of internal (closed-loop) recycling and the 
numerical example described above. The identification of 
situations in which there is a unique solution, and the way 
an allocation step affects / and/or J is the main topic of this 
paper. The paper does not discuss how the allocation prob¬ 
lem should be dealt with, although some main directions 
are indicated in the concluding section. 

This paper starts with the algorithmetic procedures of the 
inventory analysis for a descriptive LCA. Next, a number of 
very simple cases will be analyzed, with the purpose of com¬ 
ing to an exact definition of where and how computational 
problems occur that are related to the allocation problem. 
The paper closes with a summary of the findings from these 
cases, an exact definition of the allocation problem, and a 
number of implications for designing strategies to deal with 
the allocation problem. 

2 The Algorithmetic Procedures of the Inventory 
Analysis 

This paper presumes that the reader has a basic knowledge 
of the principles of life cycle inventory analysis. Neverthe¬ 
less it starts with a concise overview of selected elements of 
inventory analysis. The reason for doing so is that it concen¬ 
trates on its algorithmetic structure. Although the reader 


might feel that the algorithm for finding an inventory table 
is a dull exercise not to be discussed in methodological pa¬ 
pers, we will see that this is a crucial aspect of appreciating 
the nature of the allocation problem. Moreover, in doing so 
we are able to discover connections to other branches of 
science, in particular to economic activity analysis (see e.g. 
Koopmans 1951b), a field of science that goes back to the 
early 1950s and that has acquired a much more solid foun¬ 
dation than life cycle inventory analysis has at present. 

We assume that the economic unit processes which build 
the product system have been recorded in some format so 
that there is a clear separation between inputs from other 
processes ("economic inputs"), outputs to other processes 
(“economic outputs"), inputs from the environment 
("environmental inputs"), and outputs to the environment 
("environmental outputs"). (ISO labels economic flows as 
"intermediate product flows" and environmental flows as 
"elementary flows".) In order to construct an algorithmetic 
procedure we must first assign mathematical symbols to these 
flows. For that purpose we adopt the following rules: 

• the /th economic flow of process j is given the symbol 

• the £th environmental flow of process /' is given the sym¬ 
bol b h ; 

• input flows are represented by negative values of the vari¬ 
ables a and b; 

• output flows are represented by positive values of the 
variables a and b ; 

• flows that are not involved in a certain process have val¬ 
ues 0 for the variables a and b; 

• the number of economic flows is denoted by 1, the number 
of processes by J, and the number of environmental com¬ 
modities by K; these numbers are not necessarily equal 
to one of the others, neither is there a fixed ordering of 
these numbers in the sense that, say, K < J < /. 

Figure 1 illustrates these conventions. 


economic inflows —> 

(a v , a 2( ,.... a f .... a) —> economic 

environmental inflows -* process j 
(b, r b 2 , .V- b K )' 


-> economic outflows 
(a, ; , s 2l . a r . a)' 

environmental outflows 
-» (by, b 2i .V- bj 


Only when negative 
T Only when positive 


Fig. 1: Representation of an economic unit process 
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The above definitions comply with previous papers on the 
aigorithmetic structure of life cycle inventory analysis, in 
particular with MOllf.r (1992), Heijungs (1994), Frisch- 
knecht &C Koi.m (1995), and Heijungs (1997), although ir¬ 
relevant deviations of the chosen symbols may occur. Re¬ 
stricted to economic flews, the same structure appears in 
the much older literature of economic activity analysis (see 
for instance Von Neumann, 1945/1946; Koopmans, 1951a; 
Georgescu-Roegf.n, 1951; Chenery & Clark, 1959;Debreu, 
1959; Rosenbluth, 1968; Saxton & Ayres, 1971 and 
Ginsburgh & Waelbroeck, 1981). 

It should be observed that every process can in principle be 
multi-input and multi-output, and that representation of 
multi-function processes is not in any sense problematic. 
There is therefore no allocation problem at this stage of data 
collection, translation into an algorithm, and representation 
in a product system. 

Usually a certain functional unit has been defined in the goal 
definition of the Life Cycle Assessment. Let us assume that 
the functional unit is 1000 kWh electricity (In ISO terminol¬ 
ogy: the functional unit is electricity and the reference flow 
is 1000 kWh). The last step of the inventory analysis is to 
scale the processes that build the product system in an ap¬ 
propriate way such that the system as a whole produces ex¬ 
actly the functional unit, and that all other economic flows 
have been traced back to environmental flows. In other 
words, there is one economic output: 1000 kWh electricity, 
and there are no other economic outputs or inputs. The trans¬ 
lation of this statement into a mathematical expression which 
can be subject to an algorithm is the following: find multi¬ 
plication factors for every process such that the net economic 
flows is restricted to an output of 1000 kWh electricity, while 
all other net economic flows are zero. Denoting the multi¬ 
plication factor for process / by t, we seek to determine the 
ts such that the aggregated amount of economic flow i over 
all scaled processes equals either the magnitude of the func¬ 
tional unit for the flow that corresponds to the nature of the 
functional unit, or that it equals 0 for all other economic 
flows. If we denote the net amount of economic flow /' to or 
from the system as a (which is again negative for a flow 
entering the system and positive for a flow leaving the sys¬ 
tem), we thus have that one of these a is equal to +1000 
kWh electricity, and the other as are equal to 0. The math¬ 
ematical procedure to find the scaling factors t. for every 
process j can now be stated as: 

V i - 1. I- E Vj - «, ( 1) 

Readers that are familiar with matrix algebra will immedi¬ 
ately observe that the quantities in this equation can be sum¬ 
marized into a matrix equation: 

A • t = a (2) 

where the columns of the matrix A represent data with re¬ 
spect to one unit process, the rows of the matrix A and the 


rows of the vector a represent data with respect to economic 
flows of one single type, and the rows of the vector t repre¬ 
sent the scaling factors for the different processes. The ma¬ 
trices are thus defined as follows: 



/ 

a n 

°12 

... a}j ... 

\ 

a u 


/ \ 

“i 


/ \ 
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a ZJ 

a u 


«2 


h 

A - 

a i\ 
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... a b ... 
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The matrix A is in the economic literature often called the 
technology matrix (Van Rijckeghem, 1967; Ten Raa, 1995). 
The scaling factors t are variously called intensities (Von 
Neumann, 1945/1946), activity levels (Chenery &c Ci.ark, 
1995; Saxton &c Ayres, 1975) or operating times (Heijungs, 
1997). The set of scaling factors is sometimes called a pro¬ 
gram (Chenery &c Clark, 1959). 

Finding the scaling factors for the processes is of course not 
the aim of the inventory analysis: it is an intermediate step 
in the computation of the total environmental flows. To be 
more precise, the scaling factor of a process scales all data 
of that process, not only the data on economic flows. The 
scaled environmental flows are therefore given by the prod¬ 
uct of the recorded flow b t/ and the scaling factor for that 
process f. Aggregation of this amount over all processes yields 
the aggregated amount of environmental flow k. Denoting 
the aggregated amount of environmental flow k by f^,, we 
have symbolically 

V*. 1, .... K: p t - £V; W 

j-i 

Again we may easily identify this with matrix operations by 
putting the elements b kj into a matrix B and the elements P x . 
into a vector (3: 

P=B ■ t (5) 

Here, the matrix B and the vector p are defined as 




bn 

~ b v • 
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• b u 


/ \ 
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The vector P which is calculated here contains K environ¬ 
mental flows, often called the environmental interventions 
(ISO does not reserve a seperate term for this, and employs 
"elementary flows"). It can therefore be said to represent 
the inventory table associated with the functional unit de¬ 
fined by the only non-zero element of the vector a. Matrix 
B has, as far as we know, not been denoted otherwise than 
the intervention matrix (Heijungs, 1997). 

Once we know the scaling factors (the vector t) the calcula¬ 
tion of the inventory table (the vector p) is straightforward. 
The problem for the moment is the computation of the scal¬ 
ing factors from the process data (the technology matrix A) 
and the functional unit (the vector a). Equation (2) is an 
equation which states the relationship between A, t, and a, 
but does not express how t is to be found for a given A and 
a. Elementary matrix algebra states that t can be found by 
inversion of the matrix A and postmultiplication with the 
vector a: 

t = A 1 • a (7) 

provided A is invertible. 

It is exactly this condition of invertibility which determines 
if there is a computational problem that is very closely re¬ 
lated to the allocation problem (cf. Heijungs, 1994). If A is 
invertible, we can simply calculate t according to Equation 
(7) and proceed to calculate the inventory table according 
to Equation (5). Or in one step, substituting Equation (7) 
into Equation (5), we obtain the fundamental equation of 
life cycle inventory analysis: 


3 Some Illustrative Case Studies 

This section is devoted to a discussion of a number of cases 
in which there are certain problems of finding the inverse of 
the technology matrix and thereby finding the inventory ta¬ 
ble. It should be emphasized that the cases are completely 
fictitious, the figures are chosen such that transparent ma¬ 
nipulation is possible, instead of representing real process 
data. The first case will serve to illustrate the general for¬ 
malism that is described above. The other four cases will 
throughout be based on small modifications to the original 
process data of the first one. For easy reference, the modifi¬ 
cations will be marked by shading. 


3.1 Case i: a simple illustration without complications 

Suppose that we have two unit processes: fuel production 
and electricity production (—> Fig. 2), and that hypothetical 
process data are given (—» Table 1). The mathematical sym¬ 
bols to be assigned to these data are now as follows: a n = 
10, a n = -1, b u = 1, etc. Although we have not yet defined 
a functional unit, we may draw a system boundary around 
these two processes and thereby construct a product sys¬ 
tem. The output that is to be delivered can then be imposed 
by requiring a certain functional unit to be produced. The 
technology matrix is given by the first two rows of Table 1: 


' 10 -500] 

,-i iooJ 


(9) 


/J = B • A 1 • a 


( 8 ) 


The intervention matrix is given by the last three rows: 


If however the technology matrix A is not invertible we must 
either give up the idea of calculating an inventory table or 
we must go into a trick to find some kind of solution. 


B - 


' 1 
0.1 
, 0 


10 ' 

2 

-50, 


( 10 ) 


fuel 

electricity 

-» electricity 

electricity -* 

fuel 

-> fuel 


production 

^C0 2 

oil -» 

Production 

->C0 2 



-> S0 2 



^S0 2 


Fig. 2: The two hypothetical unit processes that play a central role in the examples 
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Table 1: Hypothetical process data for the simple example 

of two processes and two economic flows (Case i) 


economic or environmental input or output 

process 1 

(electricity production) 

process 2 

(fuel production) 

economic flow 1 (electricity in kWh) 

10 

-500 

economic flow 2 (fuel in 1) 

-1 

100 

environmental flow 1.(C0 2 in kg) 

1 

10 

environmental flow 2 (S0 2 in kg) 

0.1 

2 

environmental flow 3 (crude oil in 1) 

0 

-50 


Assuming like above the functional unit to be 1000 kWh elec¬ 
tricity, we can build a vector of external economic flows a: 

(10001 


Calculation of the scaling factors is now simply given by 


[ 10 -50o1 ( 1000l 

(-1 100 J '{ 0 J 


\ \ 

0.2 1 | 1000 

, 0.002 0.02J { 0 , 


200 
l 2 , 


( 12 ) 


Hence the scaling factor for process 1 is 200 and for process 
2 it is 2. The inventory table is easily found as 


/ 

1 

10 

/ \ 


/ 

220 

0.1 

2 

200 


24 

, 0 

-50, 

\ 2 , 


,-100, 


In other words, the environmental interventions for the cra- 
dle-to-grave production of 1000 kWh electricity are an emis¬ 
sion of 220 kg C0 2 , an emission of 24 kg SO,, and an ex¬ 
traction of 100 1 crude oil. 

We may interpret these matrix equations in terms of a sys¬ 
tem of two equations with two variables. The variables are 
the scaling factors for the two processes f, and t 2 , and the 
equations express the balancing requirement (hence the ob¬ 
solete name ecobalance) for each economic flow: 

JlOr^-500/j - 1000 (i4) 

|-l/ r 100/ 2 » 0 

The matrix formalism enables to easily find the values for 
the variables 1, and t r 

From this example two things are immediately clear: 

• The formalism has no problems in dealing with loops, 
such as that electricity production needs fuel and fuel 
production needs electricity (cf. Fava et al., 1991 and 
Heijungs et ah, 1992). 

• The same inverse of the technology matrix can be used 
to calculate any desired functional unit. As finding the 


inverse of a matrix is the most time-consuming step of 
the calculations one might consider to provide with every 
database the inverse of the technology matrix (cf. 
Frischknecht et ah, 1994). 

The next example will introduce a very small complication 
related to a process which is not involved in the product 
system providing the functional unit. 


3.2 Case ii: a process that is not involved in the product 
system 

Suppose that our database contains an additional process: 
waste incineration. The process data are given by Table 2. 
The technology matrix is now given by the first three rows 
of Table 2, and it contains an additional column for the 
newly introduced process; furthermore the intervention 
matrix also contains an extra column, and the vector of ex¬ 
ternal economic flows an extra row: 
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10 

-500 
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-5 
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10 

900 


( \ 

1000 

A - 

-1 

100 

0 

; b - 

0.1 
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10 
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-1000, 



-50 
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, 0 , 


Calculation of the scaling factors is still simple: 



10 -500 -5 

-1 
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1000 


( \ 
200 
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We see that the scaling factors for processes 1 and 2 remain 
unaffected, and that the scaling factor for process 3 is 0, 
which means that this process is not involved in the product 
system that provides the functional unit. The inventory ta¬ 
ble is therefore the same as before: 
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Tabic 2: Hypothetical process data for the simple example of three processes and three economic flows, one process and one flow being not 
involved in providing the functional unit (Case ii) 


process 1 

(electricity production) 

economic or environmental input or output 


process 2 
(fuel production) 

process 3 
(waste incineration) 

economic flow 1 (electricity in kWh) 

10 




-500 


-5 

economic flow 2 (fuel in 1) 

-1 




100 


0 

economic flow 3 (waste in kg) 

0 




0 


-1000 

environmental flow 1 (CO., in kg) 

1 




10 


900 

environmental flow 2 (SO, in kg) 

0.1 




2 


10 

environmental flow 3 (crude oil in 1) 

0 




-50 


0 

From this example we learn one additional thing: 


' 10 

-500 

-5 ] 


' iooo'' 


• The system may contain processes which are not involved 

A - 

-1 

100 

0 


0 

(18) 

in a certain product system. In other words, it is nor nec- 


0 

0 

-1000 


0 

essary to delete or "hide" processes from the database in 


0 

0 

-200 , 


0 ) 



order to calculate the inventory table. With respect to 
memory requirements of computer equipment it may be 
advisable to do so, but in the light of the previous con¬ 
clusion concerning the universal applicability of the in¬ 
verse of rhe full technology matrix, it could still be ad¬ 
vised to include all processes. 


This matrix can not be inverted because it is not square: 
there are more rows than columns. We may interpret this as 
the problem of an overdetermined system: there are three 
variables (the scaling factors for the three processes) and 
there are four equations (one for each economic flow): 


The same should obviously hold true for multi-function proc¬ 
esses that are not involved in the product system. That situ¬ 
ation, however, introduces some important complications, 
and therefore deserves a separate discussion. 


10/ r -500/ 2 .-5r 3 - 1000 
-l/,.100/ 2 *0r 3 - 0 
Of,.0/ 2 .-1000/ 3 . 0 
0/,.0t 2 .-200r 3 - 0 


(19) 


3.3 Case iii: a two-function process that is not involved in 
the product system 

Suppose that the process of waste incineration of Case ii 
fulfils two functions: processing of organic waste and process¬ 
ing of chemical waste (—» Table 3 for the process data). The 
technology matrix and the vector of economic flows now 
have four rows: 


In general, one cannot solve an overdetermined system of 
equations, because there is a conflict between two or more 
equations. 

An intuitive argument, however, says that the situation is 
not different from that of Case ii: process 3 is not active and 
columns 3 and rows 3 and 4 might have been discarded 
without altering the characteristics of the system as involved 


Table 3: Hypothetical process data for the simple example of three processes and three economic flows, one multi-function process and 
two flows being not involved in the product system that provides the functional unit (Case iii) 


economic or environmental input or output 

process 1 

(electricity production) 

process 2 
(fuel production) 

process 3 

(waste incineration) 

economic flow 1 (electricity in kWh) 

10 

-500 

-5 

economic flow 2 (fuel in 1) 

-1 

100 

0 

economic flow 3 (organic waste in kg) 

0 

0 

-1000 

economic flow 4 (chemical waste in kg) 

0 

0 

-200 

environmental flow 1 (C0 2 in kg) 

1 

10 

1000 

environmental flow 2 (S0 2 in kg) 

0.1 

2 

30 

environmental flow 3 (crude oil in 1) 

0 

-50 

0 
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for the functional unit. In fact the third and fourth equation 
are redundant: it follows directly from 

Oq + Otj+-1 OOOtj ■ 0 {20} 

0r,-0L,*-200/ 3 - 0 

that t i = 0. An important question is now how this intuitive 
aspect can be formalized so that software for life cycle in¬ 
ventory analysis is able to deal with this situation. 


gated by computing the distance between what is actually 
delivered externally (A t = A A f -a) and what was supposed 
to be delivered externally (a): 


(A-A^-a - a|| 


10 

-500 

-5 

-1 

100 

0 

0 

0 

-1000 

0 

0 

-200, 


200 

2 

0 


1 1000 N 
0 
0 
0 


. 0 


(24) 


For that purpose we need to go briefly into the theory of the 
pseudo-inverse of matrices (or Moore-Penrose generalized 
inverse; see Golub &c Van Loan, 1993 and Stewart &c Sun, 
1990 for a theoretical account, and Press et al., 1992) for a 
"numerical recipe" for calculating a pseudo-inverse. Further¬ 
more, many mathematical packages (like Mathematica, 
Matlab and Linpack) contain a standard routine for it. The 
pseudo-inverse of a matrix A is equal to the inverse for a 
normal non-singular square matrix, and is defined as well 
for rectangular matrices and for square matrices for which 
the normal inverse does not exist. The pseudo-inverse can 
be found by means of singular value decomposition, and 
has the property that it is an approximation to what a nor¬ 
mal inverse should be. While a normal inverse A' 1 has the 
property that it gives the unit matrix I when multiplied with 
the matrix itself (A A* 1 = I), the pseudo-inverse A f is the ma¬ 
trix which gives, when multiplied with the matrix itself, a 
result that is closest to the unit matrix. This should be un¬ 
derstood in terms of the systems of equations a = A t in 
which a and A are known and the unknown t can be found 
by t = A ' A t = A ' a. For a normal inverse A' 1 we have thus 


In this case the solution found with the pseudo-inverse is 
therefore an exact solution to the problem that was to be 
solved. The inventory table can be found in exactly the same 
way as was discussed in Case ii. 

One further remark here is that the system could be solved 
for this particular functional unit. The approach with the 
pseudo-inverse will for some choices for the functional unit 
not lead to an exact answer. For example, for the functional 
unit "treatment of 1000 kg chemical waste", we have 


a 


0 

0 


(,-lOOOj 


(25) 


but the solution found with the pseudo-inverse appears to 
produce an external economic flow that is given by 


H A • A' 1 • a - a II = 0 

while for a pseudo-inverse A f we have 


II A • A f • a - a II = minimal 


( 21 ) 


( 22 ) 


A-A t -a 


0 ' 
0 

-192 

,-38.5, 


The vector A'A f a - a is called the minimum residual (Golub 

& Van Loan, 1993). Its length II A A* a - a II indicates the w ' r ^ a minimum residual distance 
distance in terms of the square root of the sum of the squares 

of the different elements of the expression inside. II A ■ A + • a - a II = 196 


(26) 


(27) 


In the example above we can compute the pseudo-inverse of 
the technology matrix and postmultiply it with a to obtain 
the scaling factors: 


' 10 

-500 

-5 ' 

t 

' 1000' 


( 







200 

-1 

100 

0 


0 


2 

0 

0 

-1000 


0 



> o 

0 

-200, 


> 0 , 


, 0 , 


(23) 


Clearly we find that process 3 is not involved for the func¬ 
tional unit. An important check is now whether the pseudo¬ 
inverse has yielded the closest but yet inexact solution or 
whether it has found an exact solution. This can be investi- 


as the apparently closest solution. In the above calculations, 
we have neglected to specify the dimensions of the elements. 
A situation with non-zero minimum residuals will most of¬ 
ten exhibit bizarre unphysical units; see also Case iv. 

One new aspect has been learned: 

• Systems with multi-function processes can, when these 
processes are not involved in the product system that 
provides the specified functional unit, be solved in an 
exact way without entering an allocation procedure. 

The next case discusses the less trivial situation that a multi¬ 
function process is involved in the product system. 
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Table 4: Hypothetical process data for the simple example of two processes and three economic flows, one economic process being a multi¬ 
function process (Case iv) 

economic or environmental input or output 

process 1 

(electricity production) 

process 2 
(fuel production) 

economic flow 1 (electricity in kWh) 

10 

-500 

economic flow 2 (fuel in 1) 

-1 

100 

economic flow 3 (steam in MJ) 

1 

0 

environmental flow 1 (CO, in kg) 

1 

10 

environmental flow 2 (SO, in kg) 

0.1 

2 

environmental flow 3 (crude oil in 1) 

0 

-50 


3.4 Case iv: a two-function process that is involved in the 
product system 

Suppose that the process of electricity production coproduces 
a certain amount of steam. The process data are in Table 4. 
The technology matrix is now 


A 


10 -500 
-1 100 
, 1 0 , 


(28) 


mized has the bizarre dimension VkWh 2 +1 2 + MJ 2 . If we 
express the third flow in GJ instead of MJ, thereby replac¬ 
ing the 1 fora 3| into 0.001, we obtain quite different operat¬ 
ing times: 



(32) 


producing an external flow of 


Like in the previous case, the matrix can not be inverted 
because it is not square. We will explore the procedure with 
the pseudo-inverse. The technology matrix A can be pseudo- 
inverted and premultiplied by the vector of external flow a 
to provide the vector of scaling factors t: 


\ 

98' 
, 0 


(29) 


Calculation of the external flows which are effectively pro¬ 
duce with these scaling factors gives 


At 


980 ^ 
-98 
V 98, 


(30) 


1000 

- 0.0002 

0.2 


1 


(33) 


with a subsequently red uced minimu m residual, but a still 
unphysical dimension: VkWh 2 + r + GJ~. One might argue that 
kWh and GJ are energy units, and that 1 oil could be con¬ 
verted into energy units as well, so that everything inside the 
square root is squared energy indeed. But in a more general 
situation we will have something like ,/m /+kg 2 + m 2 +(m : ) : +hr 2 +■ • 
with MJ of energy, kg of steel, m of wire, m 2 of coating, and 
hr of labour mixed up. The point is, that minimization of 
the sum of squares of dimension-bearing quantities always 
results in a dimension-dependent minimum. 


The distance between imposed and obtained external flow 
is therefore 



( \ 

980 


( \« 

1000 

HA-A^a-al - 

-98 

- 

0 


. 98 J 


t 0 Jl 


We conclude from the above: 

• The procedure with the pseudo-inverse does not in all 
situations produce an acceptable solution. 

Only special choices of a multiple functional unit (for in¬ 
stance, 1000 kWh electricity and 200 MJ steam) resolves 
the inconsistency between the equations. 


The solution that is produced with the pseudo-inverse is 
therefore not an exact solution. Moreover, it violates the 
rules of dimension analysis, according to which different 
units may never be added. The residual that has been mini- 


The last case discusses the situation that a multi-function 
process is involved in a rather special way. 
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Table 5: Hypothetical process data for the simple example of two processes and three economic flows, one 
used flow of recycled waste (Case v) 

process producing an internally 

economic or environmental input or output 

process 1 

(electricity production) 

process 2 
(fuel production) 

economic flow 1 (electricity in kWh) 

10 

-500 

economic flow 2 (fuel in 1) 

-1 

100 

economic flow 3 (waste in kg) 

2 

-200 

environmental flow 1 (C0 2 in kg) 

1 

10 

environmental flow 2 (S0 2 in kg) 

0.1 

2 

environmental flow 3 (crude oil in 1) 

0 

-50 


3.5 Case v: a two-function process of which both functions 
are involved in the "correct" proportion 

Suppose that the electricity production process produces a 
certain amount of waste with an energetic content, and that 
the production process of fuel can use this waste. We would 
say that the fuel production process recycles waste into fuel. 
Because the waste is internally used, this situation is often 
referred to as internal or closed-loop recycling (—> Table 5 
for the process data). The technology matrix and the exter¬ 
nal flows are given by 



10 -500 


/ 

1000 

A - 

-1 100 

; a • 

0 


, 2 -200, 


t 0 , 


Again, we have a non-invertible matrix, and again intuition 
tells us that this system must have an exact solution. After 
all, the second and third equation appear to be redundant, 
because the two equations are linear combinations of one 
another: 

i0/,-500tj - 1000 

-It,. 100/, - 0 (35) 

2f,*-200t 2 - 0 

Indeed, if we employ the trick with the pseudo-inverse we 
are able to find a solution 


/ 1 


( \ 



10 -500 


1000 

( 200 

t - 

-1 100 


0 

-UJ 


i, 2 -200, 


. 0 , 


This solution can be shown to be an exact solution in the 
sense that the distance between calculated external flow and 
imposed external flow is 0. The inventory table can be cal¬ 
culated from this as usual. 


This example was based on closed-loop recycling, but the 
conclusion can be extended to systems with multi-function 
process in general. The essential element is that the func¬ 
tions that are delivered in some proportion by the multi¬ 
function process must be used in the system in exactly the 
same proportion. The situation thus also may apply to 
coproduction or combined waste treatment, not only to multi¬ 
function processes with a recycling character. Like in Case iii, 
not all choices for the functional unit will give an exact result. 
Indeed, Case iii can be seen as a special instance of Case v: the 
two functions that are delivered by the two-function process 
are involved in the "correct” proportion: 0 and 0. 

In conclusion: 

• Systems with a multi-function process can, when the func¬ 
tions of this process are involved in the "correct" propor¬ 
tion, be solved in an exact way without entering an alloca¬ 
tion procedure. Only one function is leaving the system, 
because the other functions of these processes in the sys¬ 
tem fully cover the need of the functional unit analyzed 
and therefore are entirely used up within the system. 

If a multi-function process is involved in a more general way, 
we retrieve Case iv where an allocation procedure seems to 
be necessary for finding a solution. 

This concludes our overview of hypothetical cases. It’s time 
to review the findings. 

4 Towards a Definition of the Allocation Problem 

From the hypothetical cases we are able to infer a number of 
results: 

• The matrix formalism provides a powerful algorithm for 
computing an inventory table, also for product systems 
with loops 

• The inclusion of multi-function processes leads under a 
number of different and frequently occurring circum¬ 
stances not to an allocation problem: when these proc- 
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esses turn out not to be involved in the product system, 
and when the useful flows of these processes turn out to 
be involved in the "correct" proportion 
• The matrix formalism must be extended with a routine 
to calculate the pseudo-inverse of the technology matrix 
for those cases for which the simple matrix inversion 
does not work. It may happen that the solution found 
with the pseudo-inverse turns out to be an exact solution 
of the inventory problem, especially if the multi-func¬ 
tion processes do not take part in the product system at 
hand, or if their useful flows are involved in exactly the 
"correct" proportion, such as is the case with internal or 
closed-loop recycling. 

We thus have achieved a fairly complete picture of the 
algorithmetic structure of the inventory analysis, and have 
seen a number of quite different situations in which a simple 
solution to the calculation problem of the inventory analy¬ 
sis fails. For many of those situations (such as multi-func¬ 
tion processes that are not involved and multi-function proc¬ 
esses of which the functions are involved in the "correct" 
proportion), we have found a systematic procedure. 

In the way described above, we have defined the allocation prob¬ 
lem nor in the usual way related to the occurrence of multi-func¬ 
tion processes, but to the fact that the technology matrix A can 
not be inverted, added to the fact that application of the pseudo¬ 
inverse does not yield an exact solution. So now, finally, we are 
able to define the allocation problem: 

The allocation problem refers to the fact that the scaling 
factors for the processes in the product system can not be 
determined so as to satisfy the functional unit. This condi¬ 
tion can be formalized as || A A f a - a||* 0, where A is the 
technology matrix that contains the process data, a is the 
vector of external flows that specifies the functional unit, 
A f denotes either the normal inverse or the pseudo-inverse 
of the technology matrix, and || ... || represents the square 
root of the sum of the squares of the elements (the Pythago¬ 
rean length). 

It must be emphasized that the allocation problem is hence 
not defined at a process level, but at a system level. In other 
words, one can not tell whether a particular process will 
create an allocation problem. 

Furthermore, the allocation problem is to some extent de¬ 
fined without regard to the functional unit a: the allocation 
problem is in principle present or not only in relation to the 
processes of which data is known. One might argue that the 
set of processes included in the product system depends on 
the functional unit chosen. For practical reasons that is in¬ 
deed the case. In a more fundamental respect, however, the 
set of processes included for any functional unit contains all 
processes that exist (Hofstetter, 1996). Therefore the flow 
chart of processes that build the product system is - in quali¬ 
tative terms - theoretically identical for every functional unit. 
With the development of large highly connected databases, 


like the one of Frischknbcht et al. (1994), this theoretical 
argument has already starred to become a real fact. 

In the above, it was stated that the allocation problem is 
only to some extent defined without regard to the functional 
unit. To some extent it is, however: Case iii, for instance, 
did not result in an allocation problem for the functional 
unit of 1000 kWh electricity, but did result in one for the 
functional unit of treating 1000 kg chemical waste. 

The situation is as follows: 

• If the normal inverse of the technology matrix exists, 
there is never an allocation problem (Cases i and ii) 

• If the normal inverse of the technology matrix does not 
exist because the matrix has more rows than columns, 
the pseudo-inverse may in a substantial number of cases 
provide an exact solution (Cases iii and v). There is not 
for all choices of the functional unit an exact solution, 
which indicates an allocation problem 

• If the procedure with the pseudo-inverse does not pro¬ 
duce an exact solution, there is an allocation problem 
(Case iv) 

This last point implies that it is the database itself, the col¬ 
lection of process data that is used for finding the inventory 
table for any functional unit, that can create the allocation 
problem. The allocation problem is there or is not there, 
independent of the case study at hand. However, the solu¬ 
tion of the allocation problem may be (highly) dependent 
on the goal and scope of the case study. A first attempt to 
link goals and allocation procedures is provided by 
Frischknecht (1997). 

The condition of being invertible in the normal sense can 
be translated into two requirements which both have to be 
satisfied: 

• the matrix A must be square, i.e. the number of rows 1 
must be equal to the number of columns J, i.e. the number 
of processes must be equal to the number of economic 
flows; 

• the matrix A must be non-singular, i.e. the columns of A 
must be independent vectors, i.e. no process may be con¬ 
structed by a linear combination of the other processes. 

This second requirement means in practice that one must 
take care not to load the same process twice into the data¬ 
base. If one does, there will be two columns that are exactly 
identical. The question of independence is in general diffi¬ 
cult to comprehend. In any case one should be cautious not 
to put an aggregated inventory table (say, of an electromotor) 
together with its constituent processes into the same data¬ 
base, because the aggregated inventory table is by definition 
a linear combination of the constituent processes. 

One runs into the approach with the pseudo-inverse when¬ 
ever the number of processes is smaller than the number of 
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economic flows. This can only mean that there are proc¬ 
esses that produce more than one economic flow (of course 
produce means here production of valuable commodities or 
elimination of waste). The argument is now that the occur¬ 
rence of a multi-function process is a necessary condition 
for the existence of the allocation problem. It is not a suffi¬ 
cient condition, however: Case iv and v describe a system 
with a multi-function process without an allocation prob¬ 
lem. Practitioners of life cycle inventory analysis and de¬ 
signers of databases should realize that they put too much 
effort into their work when they jump into an allocation 
procedure for every multi-function process which they see. 
It seems natural to first wait if the technology matrix is nev¬ 
ertheless invertible or if the procedure with the pseudo-in¬ 
verse gives an exact solution, and only start an allocation 
procedure for Case iv. In the situation that the normal in¬ 
verse does not exist (Cases iii, iv and v) it could be argued to 
anyhow start an allocation procedure in order to strive for a 
universally applicable database. Making the same appeal to 
universal applicability, one might argue to start the alloca¬ 
tion procedure for any multi-function process. It can be 
shown that in those cases where the pseudo-inverse leads to 
an exact solution, one would have found the same answer 
when applying an allocation procedure with an arbitrary 
allocation parameter. Practically that means that one can 
save oneself a lot of trouble by only entering the allocation 
procedure when strictly required, and not, for instance, in 
the case of closed-loop recycling. 

5 Towards a Solution of the Allocation Problem 

This paper is concerned with the nature of the allocation 
problem, not with strategies to solve the problem. However, 
understanding the true nature of the problem adds sugges¬ 
tions to the discussion how to solve it. It turns out that the 
matrix representation offers once more a suitable framework 
for analysis and discussion. Most proposals for dealing with 
the allocation problem are easily described with the matrix 
formalism. A small number of remarks will be made with 
that respect. 

The allocation problem has its origins in an impossibility to 
solve a matrix equation. We may easily identify three differ¬ 
ent strategies for solution: 

• The first strategy is to modify the question. In Case iv we 
saw such a modification: if the vector of external de¬ 
mand were to be changed from 



1000 


( \ 
1000 

a - 

0 

into a - 

0 


> 0 ; 


, 200, 


hence changing the functional unit from 1000 kWh elec¬ 
tricity to a double functional unit of 1000 kWh electric¬ 
ity and 200 MJ steam, there was no allocation problem 
any longer. If, however, we are interested in the environ¬ 


mental interventions of the different coproducts, we may 
extend this strategy with an allocation procedure to par¬ 
tition the full set of environmental interventions to each 
of the net economic flows. This strategy is used in the 
systems expansion approach to achieve an equality of 
benefits (see Fleischer 6 c Schmidt, 1996). It may be re¬ 
ferred to as "allocation at the system level". 

• The second strategy is to reduce the number of equa¬ 
tions. This amounts to reducing the number of rows of 
the technology matrix and the vector of external flows. 
In Case iv, we could merge the rows that correspond to 
"electricity" and "steam" into one row bearing the name 
"energy". Doing so will take away the redundancies in 
the set of equations. This strategy is often employed when 
a certain material is coproduced and almost the same 
material is used by another process. Or if electricity and 
steam are coproduced in, say, The Netherlands and the 
electricity is used by another process in, say, Switzerland 
in the same product system that uses the Dutch steam. The 
suggested interpretation is that "open-loop recycling is 
treated as closed-loop recycling" (see Anonymous, 1998). 

• The third strategy is to increase the number of variables. 
This amounts to introducing additional processes in the 
system. In Case iv, we might introduce an alternative proc¬ 
ess to produce steam. The inventory equation would be¬ 
come solvable, yielding a negative scaling factor for the 
newly introduced process of steam production if the func¬ 
tional unit is 1000 kWh electricity. This negative factor is 
often interpreted as "avoided production", and the whole 
procedure may be referred to as the "substitution method" 
(see Ekvai.1., 1992). An alternative form of the same strat¬ 
egy is to split the multi-function process that produces elec¬ 
tricity and steam into two independent single-function proc¬ 
esses. This would also result in a solvable inventory 
equation for a functional unit of 1000 kWh electricity, yield¬ 
ing a positive scaling factor for the electricity production 
process and a zero scaling factor for the steam production 
process. This procedure may be referred to as " allocation 
at the process level" (see Huppes, 1994). 

An interesting feature of the above categorization of strate¬ 
gies is that it differs from the description and typology that 
is currently proposed by, amongst others, ISO (Anonymous, 
1998). Our conjecture is that our proposal provides a cat¬ 
egorization that is more systematic. It shows that the con¬ 
cepts of avoided emissions is only a special case of alloca¬ 
tion of environmental interventions (both increasing the 
number of variables) and leads to a particular allocation of 
environmental interventions and requirements by referring 
to the environmental performance of competing processes 
(cf. Frischknecht, 1998a). 

The three strategies must clearly be interpreted as frame¬ 
works which may contain different elaborations. For in¬ 
stance, the third strategy contains at least two families (the 
substitution method and allocation at the process level), and 
each of these can be implemented in different ways. For the 
first family there is the choice of the most adequate process to 
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be subtracted, for the second family, the most striking differ¬ 
ence is perhaps the choice of the allocation factor on, e.g., a 
physical basis, an economic basis, marginal relationships, 50- 
50 partitioning, et cetera. Discussion of the different members 
to be distinguished within each family of strategies and of the 
ways to actually achieve the allocation procedure is, however, 
beyond the scope of this paper (see Frischknecht, 1998a, 
1998b). Suffice it to point out for the moment that strategies 
for dealing with the allocation problem go back to the eco¬ 
nomic literature on activity analysis (see, for instance. Van 
Rijckeghem, 1967; Ten Raa et al., 1984 and Konijn, 1994) 
and that a satisfactory answer has not yet been formulated 
(see Rosenbeuth, 1968 and Ten Raa, 1988). It has also been 
studied in the context of environmental economics (see 
Oenning, 1997). We think that the debate on allocation in 
LCA can be enriched by taking this tradition into account. 
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